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The classical Poisson-Boltzmann (PB) theory of electrolytes assumes a dilute solution of point 
charges with mean-field electrostatic forces. Even for very dilute solutions, however, it predicts 
absurdly large ion concentrations (exceeding close packing) for surface potentials of only a few 
tenths of a volt, which are often exceeded, e.g. in microfluidic pumps and electrochemical sensors. 
Since the 1950s, several modifications of the PB equation have been proposed to account for the 
finite size of ions in equilibrium, but in this two-part series, we consider steric effects on diffuse 
charge dynamics (in the absence of electro-osmotic flow). In this first part, we review the literature 
and analyze two simple models for the charging of a thin double layer, which must form a condensed 
layer of close-packed ions near the surface at high voltage. A surprising prediction is that the 
differential capacitance typically varies non-monotonically with the applied voltage, and thus so does 
the response time of an electrolytic system. In PB theory, the capacitance blows up exponentially 
with voltage, but steric effects actually cause it to decrease above a threshold voltage where ions 
become crowded near the surface. Other nonlinear effects in PB theory are also strongly suppressed 
by steric effects: The net salt adsorption by the double layers in response to the applied voltage is 
greatly reduced, and so is the tangential "surface conduction" in the diffuse layer, to the point that 
it can often be neglected compared to bulk conduction (small Dukhin number). 



I. INTRODUCTION 

In this two-part series, we develop a simple analyti- 
cal theory for the dynamics of electrolytes, taking into 
account steric effects of finite ion size. Motivated by 
recent experiments in microfluidics, microbatteries, and 
electrochemical sensors, our motivation is to describe the 
response of an electrolyte to an applied voltage of several 
volts, which is large enough to cause crowding of ions near 
a surface, even if the bulk solution is very dilute and in 
the absence of surface reactions. The ions thus lose their 
classical Poisson-Boltzmann distribution, which has ma- 
jor implications for their dynamics. 

As a guide to the reader, we summarize the main re- 
sults. The present "Part I" begins in this section with 
a historical review of dilute solution theory, its limita- 
tions at large voltages, and attempts to account for steric 
hindrance, specific interactions, and many-body electro- 
statics. As a first approximation, we focus only on steric 
effects, and analyze the dynamical response of a thin dif- 
fuse layer to a large applied voltage, using two simple 
continuum models (section |TT|. The key results, common 
to both steric models, are: (i) the diffuse layer's differen- 
tial capacitance is bounded and decreases at large voltage 
(section IllH) . and (it) it cannot easily engulf enough ions 
to perturb the bulk concentration or to conduct signifi- 
cant tangential currents (section IIV[) . These predictions 
are completely opposite to those of dilute solution theory 
(based on the Gouy-Chapman model of the diffuse layer). 
In the companion paper " Part 11" , we propose general, 
time-dependent equations with steric constraints and re- 
visit the parallel-plate charging problem in more detail. 



A. Dilute solution theory 

For the past century, dilute solution theory has pro- 
vided the standard model of electro-diffusion of ions [H, 
llO, mi and electrokinetic phenomena [13, [1^. The fun- 
damental assumption is that the chemical potential of a 
point-like ion i in a dilute solution has the simple form, 



fXi — kTlnCi + ZiC^ 



(1) 



where z^e is the charge, Ci the concentration, and -0 the 
electrostatic potential, determined in a mean-field ap- 
proximation by Poisson's equation. 



(2) 



typically with a constant permittivity e. The form ([T]) is 
sometimes called the "ideal" component of the chemical 
potential [37j . to which various "excess" components at 
finite concentration can be added (see below). 

In many situations, it is assumed that the ions are 
in quasi-thermal equilibrium with a Boltzmann distribu- 
tion. 



eip/kT 



(3) 



with a reference concentration c^, in which case Equa- 
tion (12) reduces to the Poisson-Boltzmann equation 
(PB). For example, to analyze a thin double layer, it is 
natural to choose c° to be the concentration of species i in 
the nearby neutral bulk solution. In most situations, the 
PB equation is hard enough to solve that the the Debye- 
Hiickel linearization for small potentials, \'4>\ <IC Vt, is 
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required for analytical progress, where tpx = kT/ziC is 
the thermal voltage (« 25 mV for monovalent ions at 
room temperature). 

The well known exception, which admits a simple so- 
lution, is the case of a symmetric binary z : z electrolyte 
in a ID geometry, where the PB equation in the form 



= a^=2zeco smh( — 



(4) 



was solved analytically by Gouy [4l| and Chapman [22| 
(GC) for a semi-infinite electrolyte of bulk concentration 



Co 



— ^0 



c_ near a flat charged surface 0,111, [t^. (It 
is less well known that Gouy [4l| also solved the case of 
an asymmetric 2z : z electrolyte.) Although this solu- 
tion may seem quite special, it describes the most com- 
mon situation where diffuse charge is confined to a thin 
capacitor-like "double layer" near a solid surface. The 
width of the diffuse layer is the Debye [IB] (or, more 
properly, Gouy [4lj) screening length. 
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ekT 
2z^e'^CQ 



(5) 



as can be seen by scaling (jH) with x = x/X]j and 
tjj = ijj/ipT to obtain the dimensionless form, tJj" — sinh ip. 
The screening length {Xd ~ 1 — 100 nm in aqueous 
electrolytes) is typically much smaller than any geo- 
metrical length, so bulk solution remains quasi-neutral 
with diffuse charge confined to thin, locally flat, quasi- 
equilibrium double layers. 

Due to its analytical simplicity and wide applicabil- 
ity, the Gouy-Chapman solution has become the stan- 
dard model for highly charged double layers. It is the 
basis for the classical theory of tangential "surface con- 
duction" through the diffuse layer [H, [ll], [g^, [13], as 
well as the recently noted phenomenon of salt "adsorp- 
tion" (or "uptake") from the neutral bulk by the dif- 
fuse layer in response to a large applied voltage [1, H^l . 
Such predictions of the GC model have major implica- 
tions for electrokinetic phenomena of the first kind, such 
as electro-osmosis, electrophoresis, streaming potential, 
and diffusiophoresis, at large surface potentials (i?!. [60|. 

Dilute solution theory has also been used in nearly ev- 
ery model of diffuse charge during the passage of current. 
Near equilibrium, the flux density of ion i is proportional 
to the gradient of its chemical potential ([T]), 



-A Vc 



kT 



(6) 



where Einstein's relation, Di — bi/kT, relates the ion's 
mobility hi to its diffusivity Di. For a system in quasi- 
steady state, V ■ Fi = 0, the nonzero current, J = 
'Y^^ZicFi, only slightly perturbs the Boltzmann distri- 
bution of the ions. As a result, the GC solution is also 
widely used to describe diffuse-layer effects on electrode 
reaction kinetics [1, It^] (the Frumkin correction [s^]). 
up to Nernst's limiting current, where the bulk concen- 
tration of an electro-active species vanishes. At larger 



"superlimiting" currents, the PB equation loses validity, 
but dilute-solution theory is still used to describe diffuse 
charge, which loses its Boltzmann distribution [9l| and 
extends into the bulk as "space charge" [s?! , while retain- 
ing an inner boundary layer at the screening length [2^. 
The dilute-solution space-charge model is the basis for 
theories of electrokinetic phenomena of the second kind, 
such as super- fast electrophoresis 29] and hydrodynamic 
instability at a limiting current [83| . 

Dilute solution theory has also been used to describe 
the dynamics of electrolytes, subject to time-dependent 
applied voltages. The classical description comes from 
the Poisson-Nernst-Planck equations (PNP), which con- 
sist of Eq. ^ and mass conservation laws for the ions. 



dt ' 



D, 



(7) 



where Di ~ constant is normally assumed. Again, 
general solutions are only possible upon linearization, 
where the potential satsifies the Debye-Falkenhagen 
equation [26|. In the nonlinear regime, for thin dou- 
ble layers, the GC solution again holds since the diffuse- 
charge remains in quasi-equilibrium, and the diffuse-layer 
acts like a voltage-dependent differential capacitance in 
series with a bulk resistor [g^]. Such equivalent circuit 
models can be derived systematically from the PNP equa- 
tions by aymptotic analysis, which also reveals correc- 
tions at large voltages, due to bulk diffusion in response 
to salt adsorption Q and surface conduction at large 
applied voltages. Such nonlinear diffuse-layer effects 
could be important in intepretting impedance spectra 
in electrochemical sensing or in understanding 

high-rate thin-film rechargeable batteries [H, HI, [H, |9^ . 

Another current motivation to study nonlinear diffuse- 
charge dynamics comes from "induced-charge" electroki- 
netic phenomena The preceding models of double- 
layer relaxation have been used extensively in theories of 
AC pumping of liquids over electrode arrays [l|, 0, SO, 
[tsI . [79I ? ], induced-charge electro-osmotic flows around 
metallic colloids [H, [66| and microstructures 0, [13, [121 , 
and dielectrophoresis ' 87l. [s ol and induced-charge 
electrophoresis 5, 25, 93, [9^| of polarizable particles in 
electrolytes, although a nonlinear analysis based on the 
PNP equations has not yet been attempted for any of 
these situations. This may be a good thing, however, 
since we will show that dilute solution theory generally 
breaks down in the regime of experimental interest, where 
the applied voltage is much larger than the thermal volt- 
age, V > IpT- 



B. Validity of the nonlinear model for dilute 
solutions 

Dilute solution theory provides a natural starting point 
to understand nonlinear effects in electrolytes, and the 
GC model is used in all of the examples above to model 
the diffuse layer at large applied potentials. In spite of its 
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mathematical appeal, however, nonlinear dilute solution 
theory has limited applicability, due to the exponential 
sensitivity of the counter-ion concentration to voltage. 
For example, in the typical case of quasi-equilibrium, the 
Boltzmann factor ([3]) brings nonlinearity to the PB equa- 
tion Q when ip > ipT, but this dependence is so strong 
that the surface potential cannot get much larger without 
invalidating the assumption of a dilute solution. 

To see this, note that there must be a maximum con- 
centration, Cmax = of counterions packed with typi- 
cal spacing a near a highly charged surface. This effective 
ion size is clearly no smaller than the ionic radius (typ- 
ically « lA), although it could be considerably larger 
(several nm), taking into account hydration effects and 
ion-ion correlations, especially in a large electric field. 
For the sake of argument, consider a very dilute bulk 
solution, Co — lO^^M, with a = 3A and z — 1. The 
nonlinear regime begins at a diffuse-layer potential 
« kT/ ze = 25 mV, but the steric limit is reached in di- 
lute solution theory at ~ IS/cT/ze = 330 mV, where 



kT kT 

= In(a^co) = In 

ze ze 



Co 



(8) 



Since the solution ceases to be "dilute" well below the 
steric limit, there is only a narrow window of surface 
potentials (~ 25 — 200 mV), sometimes called "weakly 
nonlinear regime" [6], where nonlinearity arises and the 
fundamental assumption ([T]) remains valid. 

Unfortunately, the most interesting predictions of the 
non-linear PB theory tend to be in the "strongly non- 
linear regime" , where the dilute approximation fails dra- 
matically. For example, the Dukhin number, which con- 
trols the relative importance of tangential conductivity 
in a thin diffuse layer, Cs, compared to bulk conduc- 
tivity, cr, at a geometrical scale L, has the following 
form [lS[lIl,[23, 



Du = — = 4 — smh 
aL L 



ze^D 
ikT 



(9) 



assuming the Gouy-Chapman model. This is also the 
dimensionless ratio, Ts/cqL, of the excess surface con- 
centration of ions, Fs, relative to the bulk concentration, 
so it also governs the (positive) adsorption of neutral salt 
from the bulk in response to an applied voltage [1, [6l[ . 
The gene ral derivation of Eq. ([9]) assumes a thin double 
layer [24| , but in that case [Xjj <^ L) a large Dukhin num- 
ber corresponds to situations where the steric constraint 
is significantly violated, '^d > rendering Equation 
(9) inappropriate. Similar concerns apply to other non- 
linear effects in Gouy-Chapman theory, such as the rapid 
increase of the differential capacitance (defined below) 
with surface potential. 



(10) 



which would have important implications for electro- 
chemical relaxation around conductors 0, HJ, [s^ and 
for AG electro-osmosis [1, . 



C. Beyond the Poisson-Boltzmann picture 

We are certainly not the first to recognize the limi- 
tations of dilute solution theory. Historically, concerns 
about the unbounded capacitance of a thin diffuse layer 
in the GG solution PH)) first motivated Stern [qJ] to hy- 
pothesize that there must also be compact layer of ad- 
sorbed ions on the surface, as originally envisioned by 
Helmholtz . The Stern layer is where electrochemical 
reactions, such as ion dissociation (setting the equilib- 
rium charge on a dielectric) and/or redox couples (set- 
ting the Faradaic current at an electrode), are believed 
to occur, within a molecular distance of the solid sur- 
face [5^ Is^l- The Stern layer capacitance (see below) 
helps to relieve the overcharging of the diffuse layer in 
Gouy-Ghapman theory, but, due to steric constraints, it 
too cannot conceivably withstand a voltage much larger 
than Instead, at larger voltages, the region of ion ac- 
cumulation must inevitably extend away from the surface 
into the solution, where ions undergo hindered transport 
in a concentrated solution without having specific inter- 
actions with the solid. 

The most basic aspect of a concentrated solution is 
the finite molecular length scale, gp. Over half a cen- 
tury ago, Wicke and Eigen [3l|, [S^, [93| made perhaps the 
first attempts to extend dilute solution theory to account 
excluded volume effects in a simple statistical mechani- 
cal treatment. The theory was developed further in the 
past decade by Iglic and Kral-Igli c [ll. [isl jisl . Isi ] and 
Borukhov, Andelman, and Orland [l7, 1^19l|. These au- 
thors proposed free energy functionals, based on contin- 
uum (mean-field) approximations of the entropy of equal- 
sized ions and solvent molecules, which they minimized 
to derive modified Poisson-Boltzmann equations (MPB). 
The motivation for this work was mainly to address the 
effect of large ions, whose sizes may be comparable to 
the screening length, in an equilibrium diffuse layer with 
^'d ~ V'T- Our main point here is that crowding effects 
can also be very important for small ions near a polar- 
izable surface, when subjected to a "large" voltage the 
exceeding the threshold ^'c. No matter how dilute is the 
bulk solution, a sufficiently large electric field can always 
draw enough ions to the surface to create a concentrated 
solution in the diffuse layer. 

There have also been attempts to go beyond the mean- 
field approximation of steric effects, by treating specific 
(and in some cases also Goulombic) ion-ion and ion-wall 
interactions. Simple MPB equations have been proposed, 
which modify the charge density to account for the spa- 
tial correlation function between an ion in solution and 
a fiat wall (via an effective external potential) based on 
molecular dynamics simulations (33l . [sol . [76j . but such 
models are not easily extended to any other geometry, 
such as a rough wall [5^. Gorrections of this type can 
also be obtained from a general probabilistic model of 
interacting ions, whose dynamics are given by nonlin- 
early coupled Langevin equations. Using this approach, 
Schuss, Nadler and Eisenberg rigorously derived "condi- 
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tional PNP" (and PB) equations where each ion concen- 
tration Ci(r) in the mean-field Poisson's equation ^ is 
replaced by the conditional probability density of find- 
ing an ion at a certain position, given the positions of 
the other ions [13, [H, Ha] , although a simple closure of 
the model requires further assumptions about statistical 
correlations. 

There are a variety of general statistical mechanical ap- 
proaches from liquid state theory which have been 
applied to electrolytes, taking into account not only steric 
effects, but also many-body electrostatic correlations (see 
below). Since the 1970s, the modest, but challenging, 
goal has been to accurately predict the equilibrium dis- 
tribution of ions in Monte Carlo simulations of the "prim- 
itive model" of charged hard spheres in a homogeneous 
dielectric continuum bounded by a hard, charged wall. 
Typically, the model is studied in the limit of "small" sur- 
face potentials (^'^ « kT/e) relevant for equilibrium sur- 
faces. For example, a plethora of MPB equations (such as 
"MPB4", "MPB5",...) perturbing Gouy-Chapman the- 
ory have been derived by variations on the Mean Spher- 
ical Approximation [2l|, [H, [H, III [ll III Hi- More 
complicated, but often more accurate, theories have also 
been derived using statistical Density Functional Theory 
(DFT) [38,, 39, 49]. By writing the ion flux as the gradient 
of an electrochemical potential obtained from the DFT 
free energy functional [371 do in Part II, using the 

much simpler functional of Borukhov et al. [l8|), it has 
been shown that the selectivity of ion channels can be 
understood to a large extent in terms of entropic effects 
of charged hard spheres [13] . 

In spite of many successes and the appealing first- 
principles nature of liquid-state theories, however, they 
are not ideally suited for our purpose here of capturing 
time-dependent, large- voltage effects in a simple model. 
Liquid-state theories are not analytically tractable and 
thus not easily coupled to macroscopic continuum mod- 
els, although perhaps empirical fits could be developed 
for quantities of interest. For example, the DFT ap- 
proach requires nontrivial numerical methods, just to 
evaluate the steady flux in a one-dimensional ion chan- 
nel [13, [13] ■ Another concern is that liquid-state theories 
have been designed and tested only for "small" voltages 
{^D ~ kT/e) at equilibrium surfaces, so it is not clear 
how well they would perform at large, time-dependent 
voltages, since there are no rigorous error bounds. For ex- 
ample, the double-layer differential capacitance in modi- 
fied PB theories often increases with voltage [zl] , some- 
times even faster than in PB theory [13], but we will 
argue below that it must decrease at large voltages, es- 
pecially in concentrated solutions. In light of the scaling 
e'^u/kT from PB theory ([3]), a related problem is that, 
until recently (49i . i55.] . most theories have been unable 
to predict the decay of capacitance at low temperature 
in concentrated solutions. Although liquid-state theories 
for large voltages should certainly be developed, we will 
focus on much simpler mean-field, continuum descrip- 
tions of steric effects, in the hope of at least capturing 



some qualitative features of nonlinear dynamics analyti- 
cally, across a large range of applied voltages 

Consistent with this approach, we will also work with 
the mean-field continuum description of electrostatic in- 
teractions (12), which neglects discrete many-body corre- 
lations. In passing, we point the reader to some of the 
extensive literature on corrections to PB electrostatics, 
for the simplest possible model of an equilibrium liquid of 
point-like charges in a uniform dielectric medium. In the 
absence of specific interactions (such as steric repulsion) , 
the fundamental length scale for ion-ion correlations is 
the Bjerrum length (w 7 A in water at room tempera- 
ture). 



inekT 



at which the bare Coulomb energy balances the ther- 
mal energy. Interesting consequences of many-body elec- 
trostatics, not present in the mean-field PB descrip- 
tion, include Oosawa-Manning condensation of counteri- 
ons around charged rods in polyelectrolytes [^ , [^ [tI] , 
renormalized charge of colloidal spheres [1,M1 > W' 
hanced counterion localization around discrete surface 
charges , and counterion- mediated attraction of like- 
charged objects [H m, [13] • The latter phenomenon is 
believed to be responsible for the condensation of DNA 
in multivalent electrolytes [l^] , as well as the adhesion of 
cement plaste [13, [E^j . A key part of the physics is the 
attraction between an ion and its "correlation hole" re- 
sulting from a fluctuation, which has recently been incor- 
porated into a modified PB equation for a flat wall [ssl ]. 
In all of these problems, however, the equilibrium surface 
charge is typically small (up to a monolayer of ions); it 
would be interesting to study electrostatic correlations at 
a much more highly charged surface, such as an electrode 
applying a large voltage, V ^ ^pT (our focus here). 

Finally, we mention solvent effects, which are much 
less studied, and surely also very important at large volt- 
ages. Electrochemists widely believe that water dipoles 
in the Stern layer are so highly aligned by large electric 
fields that the effective permittivity drops by an order 
of magnitude (e.g. from SOeq to 5eo) [ij]. At large ap- 
plied voltages, where typical fields are of order V/nm, 
it is reasonable to expect that the reduced permittivity 
would extend into the diffuse layer. This could have a 
major effect on ion-ion correlations, since the Bjerrum 
length If, cx e^^ could get as large as 10 nm. Other as- 
pects of water structure, such as hydrogen bonded net- 
works, could also be altered by large electric fields and 
large ion concentrations. It would be interesting to per- 
form ab initio quantum-mechanical simulations of highly 
charged double layers to investigate such effects, beyond 
the primitive model. 



5 



D. Scope of the present work 



In spite of the considerable literature on MPB descrip- 
tions of electrolytes in equilibrium or in steady state con- 
duction, we are not aware of any attempt to go beyond 
dilute solution theory (PNP equations) in analyzing the 
dynamics of electrolytes in response to time-dependent 
perturbations, such as AC voltages. Accordingly, here we 
develop only some very simple models with the goal of 
identifying generic new features of diffuse-charge dynam- 
ics in a concentrated solution. As such, it is preferrable 
to start with equations that capture the essential physics, 
while remaining analytically tractable, at least in some 
representative cases. For this reason, we focus on mean- 
field theories of steric effects and specifically build on 
the MPB equation of Borukhov et al., which can be inte- 
grated analytically in a few simple geometries [itI. [Tsl. [l9j . 
Such models also make reasonable predictions across the 
entire range of voltages. 



The contribution is broken into two parts. Here, in 
Part I, we consider the canonical problem of charging 
a thin double layer, viewed as an effective circuit ele- 
ment fd\ . We begin in section |TT] by describing two sim- 
ple models for steric effects in a quasi-equilibrium diffuse 
layer: (i) A "composite layer" model, consisting of a di- 
lute PB outer region and, at high voltage, an inner "con- 
densed layer" of ions packed at the steric limit, and (ii) 
the MPB model of Borukhov et al., which decribes the 
same transition in a continuous fashion. In section IIIIl 
we then analyze the diffuse-layer capacitance and its role 
in equivalent circuit approximations. In section IIVI we 
calculate steric effects on salt adsorption and tangential 
surface conduction by the diffuse layer, and discuss how 
these high- voltage effects affect and modify the applica- 
bility of circuit models. In section |Vl we also briefly 
discuss the effects of a compact dielectric layer (which 
could model a Stern layer of adsorbed ions, or, more ac- 
curately, a thin film coating the solid), before concluding 
in section IVll 



In Part II, we consider explicitly time-dependent prob- 
lems with a general formalism, not only applicable to thin 
double layers. We start with the free energy functional 
of Borukhov et al. and derive modified Poisson-Nernst- 
Planck equations (mPNP), based on a simple generaliza- 
tion of chemical potential ([T]) for concentrated solutions. 
As an illustration, we then repeat the nonlinear asymp- 
totic analysis of Ref. for the response of a blocking 
electrochemical cell (no reactions at the electrodes) to a 
suddenly applied voltage, to expose some general conse- 
quences of steric effects at large voltage. We also clarify 
the range of validity of the thin-double-layer circuit ap- 
proximations derived here in part I. 



II. TWO MODELS OF STERIC EFFECTS IN A 
THIN DIFFUSE LAYER 

We focus on the response of a thin diffuse layer to an 
applied voltage, where it suffices to consider only quasi- 
equilibrium properties @, as we justify carefully in Part 
II. Following Gouy and Chapman, we consider the case of 
a symmetric z : z electrolyte here, although our reasoning 
is readily extendible to the general case. We also assume 
that the permittivity e is constant in space, which is cer- 
tainly not correct when dense layers of ions form close 
to the surface. However, this can be taken into account 
in a following step and does not change the qualitative 
picture emerging from the following analysis. We will 
return to this point below in section IVl 

There are at least three important lengths in our mod- 
els. The first is the Debye length Ad given by (O, which 
sets the width of the diffuse layer at low voltage and low 
bulk concentration, cq. The second is the mean spac- 
ing of ions in the bulk electroylte, — (2co)~^/'^, and 
the third is the mean spacing of ions (essentially of the 

same sign) at the maximum concentration, a = Cmlx ■ 
A fourth scale L would characterize the geometry, as in 
Part II, but here we consider the regime of thin double 
layers, where Ajj <C A fifth scale would be the Bjer- 
rum length Zs, which we neglect by making the usual 
mean-field approximation. 

From the first three lengths, there are two dimension- 
less groups. The first is a/Xo, which we assume to be 
small for simple electrolytes in somewhat dilute solutions, 
so that steric effects are important only very close to the 
surface, at the inner portion of the diffuse layer (and even 
then, only at large voltages). The second dimensionless 
group can be written as the mean volume fraction of ions 
in the bulk, 

u = 2a^co = {a/lo)\ (11) 

a natural measure of non-diluteness, which controls the 
importance of steric effects, along with the dimensionless 
voltage, ze'^o/kT. In the figures below, we display re- 
sults for V = 0.00005, 0.005, and 0.5 to span the range 
from dilute to highly concentrated solutions. 

We stress that the phenomenological parameter a is 
not necessarily the diameter of an ion, ao ~ lA. We 
prefer to think of it as a cutoff for the unphysical diver- 
gences of PB theory, which we know must exist, and our 
goal is to understand its consequences. This cutoff length 
could include at least a solvation shell, e.g. a « 3A, but 
ion-ion correlations could effectively increase it further, 
since the Bjcrrum length is at least 7 A. As noted above, 
taking into account the decrease of permittivity (by as 
much as a factor of ten) or other solvent effects at large 
fields could make lb, and thus perhaps also a, as large 
as 10 nm. As a guide to using our results, we refer to 
Fig. [T] for the value of the dimensionless parameter v for 
different values of a and cq. 
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FIG. 1: The dependence of the dimensionless parameter u as 
a function of the bulk concentration cq and effective ion size 
a. 



A. A composite diffuse-layer model 

In this model, we assume that the concentration fields 
are governed by Boltzmann distributions 



c± = coe 



(12) 



wherever they are meaningful, that is, whenever these 
concentrations are weaker than a physical limit l/a"^, 
which is set by the ion size. For both ion species, if 
the formula (fT^ yields a quantity bigger than 1/a^, we 
set the concentration field of the counter- ion to be 1/a^, 
and assume that the colons are excluded from the corre- 
sponding condensed layer. The basic physics is shown in 
Fig.H 

For most geometries, such as rough surface, this model 
is ill-posed and would require additional constraints to 
determine the location the sharp free boundary separat- 
ing the "dilute" Boltzmann region and the "condensed" 
region at the maximum density. For a flat, semi- infinite 
diffuse layer with a given total charge qo or total volt- 
age ^D, however, the transition occurs at a single well 
defined position, x — Ic, away from the solid surface (at 
a; = 0). In that case, we separate the diffuse-layer into 
two parts if the potential is strong enough 

c± = coe=F^«'A/fc^ if a^c± < 1, x > h 



c± = a 



, Czp if a^coe^^^'^^'''^ >l,0<x <lc 



(13) 

Note that in this simplistic dichotomy the concentration 
field for co-ions is slightly discontinuous at the critical 
potential mentioned earlier. 



kT , 2 
— m — 

ze V 



(14) 



where again v = 2a^CQ is volume fraction occupied by 
(all) ions at zero potential in the bulk. Although very 
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FIG. 2: (a) Most prior work in electrokinetics has dealt with 
surfaces of pre-existing equilibrium charge at the scale of one 
electron per surface atom or less. This charge can be screened 
by roughly a monolayer of ions (partly in the diffuse layer) 
which corresponds to a "small" double layer voltage, of order 
il>T = kT/e. (b) In contrast, nonlinear electrokinetics deals 
with polarizable (mainly metal) surfaces, where much higher 
surface charge densities can be produced by an applied voltage 
or electric field nearby, and thus the double layer can "over- 
charge" to the point where diluite-solution theory no longer 
applies. The existence of a minimum ion spacing implies the 
formation of a condensed layer of counterions near the surface. 



simple, this model already captures steric effects at large 
voltages to a great extent. 

Let us compute the thickness Ic of the corresponding 
layer of "condensed" counter-ions. We consider an elec- 
trode to which a strong negative potential vI/d is applied 
to the diffuse layer with respect to the bulk, such that 

d\ — —"^ D > *c which leads to a condensed layer of 
positive ions in its vicinity. 

Poisson's equation for the thin diffuse layer reads 



dx^ 



(15) 



since (for now) the permittivity e is assumed constant. 
Within the condensed layer, we have p = zec^ = ze/a^, 
so by integrating, we obtain 



dx 



ze qci , 

txH , tp 

ea-^ e 



2 ea-^ e 



where —qd is the surface charge density on the electrode, 
so that qd is the total charge per unit area in the diffuse 
layer. 

Within the outer diffuse layer, the standard PB equa- 



7 



1.5 



0.5 



v = 0.5 



V = 0.005 



V = 0.00005 



40 



60 



80 



100 



ze I'PJ / kT 



FIG. 3: The thickness of the condensed layer thickness Ic as 
a function of the total voltage drop 4' d across the diffuse and 
the condensed layers. 



tion holds, 



(17) 



which gives 



d4 
dx 



= -2 



2kTcQ ( zeij) 



(18) 



At the interface between the condensed layer and the 
diffuse layer, we require the continuity of the electric field. 



and of the electrostatic potential, so that 




(19) 



kT 



ln(2/i^) 



I* 



1 ze 



D\ 



2ea 



ll + ^lc (20) 



3 c 



Combining these equations and solving for Ic yields 



V 

2 



^ kT 



- ln(2//^)] 
(21) 

which is plotted in Fig. [3] for several values of v. Gen- 
erally, the condensed layer forms when the diffuse layer 
voltage 4" D becomes only a few times the thermal voltage 
kT / ze, and then it grows sublinearly, proportionally to 
the square root of the potential drop as anticipated from 
Poisson's equation with a constant charge density. 



From this equation and ((T9)) . we finally obtain the 
charge/potential relation for the composite diffuse layer: 



ge/ = 2zecoXD^l^ (1 - + - HV'^)] 

(22) 

which holds for \'i>D\ > ^c, as assumed here. For weaker 
potentials, there is no condensed layer and the standard 
PB model holds: 

/ ze^ D \ 

qpb = -AzecoXi, sinh ( j (23) 

We note that, when compared to the PB charge density 
Qpb, the composite-layer charge density qd is significantly 
reduced by the steric effects at higher In particular, 
it increases only sublinearly, as opposed to exponentially 
in the PB model. 



B. The modified PB model 

The second model we consider is the classical mean- 
field description of steric effects in equilibrium mentioned 
in the introduction; we refer the reader to the litera- 
ture for its derivation and statistical mechanical assump- 
tions [11 [H, [13, [11 [H, H im, m, H H p. For a 
symmetric z : z electrolyte, the concentrations in the dif- 
fuse layer as a function of the electrostatic potential with 
respect to the bulk ip are given by the modified Boltz- 
mann distribution 



c± = 



coe 



^zeip/kT 



l + 2i.sinh2 (fflj 



(24) 



where the packing parameter u = 2a^CQ is again the bulk 
ion density scaled to its maximum value and a is the 
effective size of the ions and the solvent molecules. Note 
that the concentration of each ion saturates and cannot 
exceed the steric limit. 

In a mean-field approximation with these ion con- 
centrations, the potential satisfies the modified Poisson- 
Boltzmann (MPB) equation. 



zeco 



2 sinh (^) 



^ 1 + 2i/ sinh^ 



zeip \ 
2kT J 



(25) 



Unlike the composite layer model, the MPB model can be 
applied to any geometry (just like the PB model). In the 
case of a fiat diffuse layer, it gives a similar description, 
except that steric effects enter smoothly with increasing 
voltage, and there is no sharply defined condensed layer. 

As for the first model, we can integrate the MPB equa- 
tion across a thin double layer to obtain the normal elec- 
tric field at the inner part of the diffuse layer 



dx 



-sgn(V') 



2zecoA_D 2 



ln[l -I- 2i/sinh' 



/ zeijj \ , 



\2kT)^ 



(26) 
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FIG. 4: The diffuse layer charge given by PB, MPB and CDL 
models as a function of the across potential drop. 



Integrating ([T5]) using ([^5]) . we obtain for this model the 
relation between the charge per unit area in the diffuse 
layer qmpb and the potential drop across it, ^'^j: 



Qmph 



-sgn(5'£))2zecoA£) W — ln[l + 2z/sinh^ 



D 



2kT 



(27) 

This formula is illustrated in Fig. 21 and compared 
to the analogous formula for the composite diffuse- 
layer model ((22)) . Although the MPB concentrations 
have been analyzed carefully by previous au- 
thors [11 [H, [13, [II m, d, Hi, Hi, it seems the 
charge-voltage relation (|7f|) has been overlooked, or at 
least is not stated explicitly in these papers. (Similarly, 
Chapman was the first to explicitl y w rite down the 
formula for the total charge [23], even though 

Gouy had thoroughly analyzed the concentration and 
potential profiles in the "Gouy-Chapman model" a few 
years earlier [ilj.) 

We note in passing that the asymptotic behaviour at 
large voltages is similar for the two models presented, 
and corresponds to a layer of essentially constant charge 
density, which results using Poisson's equation in a total 
charge proportional to the square root of the potential 
drop. 



III. EQUIVALENT CIRCUIT ELEMENTS 
A. A word on circuit models 

As widely recognized in electrochemistry [H, [H, [qO] , 
it is often appropriate to describe the effect of the dou- 



ble layers on the dynamic properties of a system through 
effective resistances and capacitances. The underlying 
assumption, is that the equilibration time of the double 
layer is fast compared to the dynamics of the global " cir- 
cuit" considered, essentially because it is so thin so that 
transport processes within the double layer are rapid. A 
mathematical justification, starting from the simple PNP 
equations, can be given in terms of an asymptotic anal- 
ysis of the thin double layer limit 

Therefore, we anticipate that for the many situations 
where the double layer is thin compared to the system 
size, the dynamics can be understood to a significant ex- 
tent using the equilibrium characteristics of the double 
layer. These then provide appropriate boundary condi- 
tions for the dynamic transport processes in the bulk. 

In particular, the analysis presented in Ref. shows 
that the capacitance, or more precisely the differential 
capacitance "C" , of the double layer is a central quan- 
tity allowing the modeling of the system in terms of an 
"RC" circuit. The second quantity of relevance, as it 
characterizes the entrance into the strongly non-linear 
regime, is the neutral salt uptake by the double layer 
which can result in an appreciable depletion of the bulk, 
leading to modifications of its conductivity (affecting the 
resistance "i?" in the circuit) and thus of the dynamics. 
This paper and a subsequent one 0] also pointed out (in 
the context of dilute solution theory) that the tangential 
conduction through the diffuse layer is intimately tied to 
neutral salt adsorption, and indeed is governed by the 
same dimensionless "Dukhin number" . 

Therefore, we will proceed to compute all of these dy- 
namical quantities for a thin quasi-equilibrium double 
layer, focusing on the general consequences of steric ef- 
fects, which are common to the two models. After that, 
we will return to the question of surface capacitance and 
consider the effect of a thin dielectric layer (such as an 
oxide coating on a metal electrode, or perhaps a frozen 
Stern layer of adsorbed ions) on the overall dynamical 
response of the double layer. 



B. Total and differential capacitances 

The total capacitance of the diffuse layer can be ob- 
tained directly from the previous equations relating q 
(the charge per unit area in the diffuse layer) to 
(the voltage drop across the double layer), as simply 
~q{'^ d)/'^ D- We have already computed these quan- 
tities above, and they are compared to the PB result 
{ly — 0) in Fig. [S] It is immediately obvious that the 
capacitance is greatly reduced at high voltage, compared 
to the predictions of dilute solution theory. The effect is 
so dramatic that the capacitance in both models reaches 
a maximum not much larger than the zero- voltage value, 
and even decreases for all values of voltage at sufhciently 
high concentration - the opposite trend from PB theory. 

As noted above, the differential capacitance, defined 
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ze l-PJ/kT 

FIG. 5: The total capacitance —q{^D)/'i/D of the difluse layer 

as a function of '^d- FIG. 6: The diffuse layer differential capacitance as a function 

of the potential drop across itself. 



for the diffuse layer as 



dq 



D 



(28) 



is the relevant quantity for the dynamical response to 
an applied voltage. Throughout this paper, to be clear, 
we will use the notation C only for the differential ca- 
pacitance. In the PB model {v = 0), the differential 
capacitance from Eq. ((23|) has already been noted above: 



C 



D 



- — cosh — — — 
Ad V 2fcr 



(29) 



For the composite diffuse layer (CDL) model introduced 
above, dH]) yields 



1 



D 



I2v 



4^ 



kT 



H2/iy)] 



(30) 



when > Otherwise, the PB formula ([29]) still 
holds. For the MPB diffuse layer model, using (|77|) . we 
find 



fTN-h(^)| 



[1 + 2^sinh2 (f^)] + 2^sinh2 (^)] 

(31) 

The different models are compared in Fig. [6] and FiglT] . 

The qualitative trends in C are again similar in both 
models of steric effects and dramatically different from 
PB theory. For both models, the differential capacitance 
increases at very low potentials (comparable to the ther- 
mal voltage, kT/ ze) as it does in the PB model because 
steric effects are still negligible. These enter the picture 
at larger potentials and limit the storing capacity of the 
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FIG. 7: Same as figE] except with a linear scale on the y-axis. 
The case u = Q.h is not shown. 



layer, with a differential capacitance that actually de- 
creases to zero at large potentials. As a consequence both 
models predict a non-monotonic differential capacitance, 
and show a maximum at intermediate values. 

Of course, there are some clear differences in the de- 
tailed predictions of the two models, shown in Figsl6]and 
[71 Although the limiting behaviors at large and small 
voltage are similar, the transition is unphysically sudden 
in the CDL model, compared to the more reasonable, 
smooth transition in the MPB model. This is especially 
true at low concentrations, where the sudden, delayed ap- 
pearance of steric effects in the CDL model gives rise to a 
sharp cusp in the differential capacitance versus voltage. 
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At high concentrations, where the "low voltage" regime 
of dilute solution theory effectively vanishes, the CDL 
model also fails to predict the immediate onset of steric 
effects, even at small voltages, which causes a monotonic 
decrease in differential capacitance - the opposite trend 
of the PB model. Nevertheless, the CDL model manages 
to approximate the trends of the MPB model well, with 
an appealingly simple physical picture. 

In summary, three basic features show up in both mod- 
els, which we take as an indication that they qualita- 
tively hold irrespective of the specific approximations cm- 
bedded in each model: (i) the differential capacitance 
Cd(^_d) is weaker at moderate and high potentials than 
if steric effects are neglected (as in the PB scheme); (ii) 
at moderate concentrations, the differential capacitance 
varies non-monotonously with a peak at intermediate 
voltages and a slow decrease towards zero at higher volt- 
ages, (iii) at the steric limit {ly = 1), the differential ca- 
pacitance is a strictly decreasing function of voltage (in 
the MPB model). These effects are all explained by the 
strong tendency of ions to form a condensed inner layer 
at high voltage and/or high concentration, when steric 
effects are taken into account. This greatly reduces the 
differential (and total) capacitance compared to classical 
PB theory, which neglects the finite size of ions, and thus 
predicts an absurd exponential pile-up of ions extremely 
close to the surface (less than one molecular radius) in 
the nonlinear regime. 



C. Diffuse-layer charging dynamics 

These revelations have important consequences for our 
understanding of double-layer charging in many situa- 
tions. In the simplest picture of an equivalent RC circuit, 
with R the resistance of the bulk, these statements relate 
to the response time of the system to a step or an AC ap- 
plied voltage. The typical response time for a driving of 
amphtude V is Tc{V) = RCd{V), so the classical picture 
from PB theory (I29p has been that nonlinearity greatly 
slows down the charging dynamics 2A, ,881, ^9^] . 

Although this may occur in a dilute solution for a rela- 
tively small range of voltages (typically only several times 
the thermal voltage), steric effects in a concentrated so- 
lution bound the relaxation time at a value much less 
than expected from the PB model, in both the CDL and 
MPB models. This is clear in Fig. [U where we solve the 
RC circuit dynamics 



D 



dt 



y - * 



D 



R 



(32) 



for the three models numerically, to obtain the diffuse- 
layer voltage ^D{t) in response to a suddenly apphed 
voltage V across the layer in series with a bulk resistance 
R. In addition, importantly, the response time of an 
electrolytic cell is found to decrease with the amplitude of 
the applied voltage above the threshold for strong steric 
effects, V > ^E'c- As shown in Fig[5]Jb), the relaxation is 
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FIG. 8: DifFuse-layer relaxation in response to a sudden 
voltage step across a blocking cell of width 2L and con- 
stant resistance R. Time is scaled with the charging time 
Tc = \dL/D = RC{-^D = 0) (5 = 0). (a) Relaxation of 
the diffuse- layer voltage ^I^o in the PB, MPB, and CDL mod- 
els for the same applied voltage, V = IQkT / ze. (b) Relax- 
ation of the total diffuse-layer charge for different voltages, 
zeV/kT = 1, 10, 20, in the MPB model with v = 0.005, show- 
ing varying response times due to the non-monotonic voltage 
dependence of the differential capacitance (inset). 



faster for V = 2QkT/ze than for 10fcT/2;e, since vl^c = 
8.3kT/ze for v = 0.005. 
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IV. BEYOND CIRCUIT MODELS 
A. Diffuse layer salt adsorption 

As recently pointed out in Ref. [6| , circuit models can 
break down at large voltages if a large number of ions 
(mostly counter-ions) are engulfed by the diffuse layers 
with a resulting depletion of ions in the bulk. The total 
salt concentration in the diffuse layer (counter-ions plus 
co-ions) increases with voltage, regardless of the sign of 
the charge. Therefore, a diffusion layer forms and relaxes 
into the neutral bulk whenever a voltage is applied across 
a double layer at a blocking surface (although reactions 
and/or rejection of adsorbed ions from the Stern layer 
could lead to negative adsorption or salt expulsion in 
other situations |6l|). In the absence of salt injection by 
the surface, the positive adsorption of neutral salt by the 
diffuse layer is present in the PB description, where the 
counter-ion concentrations increase exponentially with 
voltage It is still present but obviously weaker in 
models accounting for steric effects, which severely limit 
the capaticity of the diffuse layer to store additional ions 
at high voltage. We now quantify this statement using 
the two simple models introduced above. 

Following notations introduced in Ref. Q, we define 
the excess neutral salt in the double layer as 



w — w 



bulk 



(c+ + c_ — 2co)dx. 



surface 



For the PB model one finds 



10"=° ScqXd sinh^ 



(33) 



For the CDL model, the same equation holds for 5'd < 
while above this value ^f^i > 'I'c we obtain 



ui" = ( - - 2 I cqIc + 2coAd 




(34) 



with Ic{^d) to be extracted from (HI]). For the MPB 
model 



w 



''T (coshu — 1) 2co\d {I - i')du 



l + 2;ysinh u . l\r,(l + 2iysmh 



(35) 



which we can compute numerically. 

Unlike the PB formula ([33]) . which predicts expo- 
nentially diverging salt uptake for increasing "^d, the 
steric modified formulae ((34l) and ([35| predict sublinear 
(square-root like) dependence on "i'n, at large voltages, 
as can be seen by inspection of p3 |l - ([35|l or from the plots 
in Fig. [HI This can be understood qualitatively as a con- 
sequence of the roughly square-root voltage dependence 
of the condensed layer width (due to its constant charge 
density), since most of the adsorbed ions are condensed 
counterions at large voltages. (See Fig. [3] and Eq. (I^T]) 
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FIG. 9: The diffuse layer neutral salt uptake w as a function 
of the potential difference d- 



for the CDL model.) Another way to see it in both mod- 
els is that the total salt adsorption is asymptotic to the 
total charge, w ^ q/ze^ which is true for any nonlinear 
double layer model (PB, MPB, CDL,...) since the colon 
density goes to zero across most of the diffuse layer at 
large voltages. As explained in Ref. @, in response to 
an applied voltage the neutral bulk electrolyte becomes 
depleted, as the electric field draws counterions into the 
diffuse layer and conducts colons away through the bulk, 
resulting in a slowly expanding diffusion layer, reducing 
the bulk concentration accordingly. 

The main feature, again common to both models, is 
that the salt uptake at large voltage is greatly reduced 
in comparison to the exponential growth predicted by 
the PB picture. This is important for the dynamics as 
this quantity sets the limit of applicability of the widely 
used RC circuit model This equivalent circuit ap- 
proximation should thus hold up to much larger values 
of potential when steric effects are included. 



B. Breakdown of circuit models 

With analytical expressions for total salt adsorption 
by the diffuse layer, we can estimate the upper limits 
on the applied voltage where the circuit approximation 
breaks down in the various models. For an electrolyte 
cell of thickness 2L, the salt uptake by the diffuse layer 
corresponds to a removal of 2w charge carriers (ions) per 
unit surface, from a bulk that contained 2cqL such carri- 
ers initially. As long as the ratio of these two quantities 
w/{coL) is small, the total resistance of the RC circuit 
is roughly unaffected by the salt adsorption. So, an esti- 
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mate for the hmit is fd\ 



as 



CqL 



« 1, 



(36) 



which translates an upper bound on the applied voltage, 
l^^l < Vthreahoid for the RC description to remain valid 
(in the thin double layer limit \d/L <^ 1). For the PB 
model, the upper bound. 



Vt 



2kBT 



threshold 



In 



L 
4A7 



(dilute) 



is not much larger than the thermal voltage, due to the 
weak logarithmic dependence on L/Xjj. 

For the models accounting for steric effects, however, 
the upper bound is greatly increased in concentrated so- 
lutions. 



reshold 



L 



o^Cq (non-dilute) 



This shows that the widely used circuit approximation 
for a thin double layer ^ L does not break down 
until enormous voltages, oc (L/Ad)^, in a concentrated 
solution. In a very dilute solution, where v = 2a^co <C 1, 
the circuit approximation may break down at moderate 
voltages, but only in a microsystem, where the double 
layers are not so thin. For example, even in an (aqueous) 
electrolyte with A^ = 10 nm, cq = 10~^ M, a = and 
in a microdevice with L = lOO^m features, the threshold 
voltage (with steric effects) is roughly 0.2 Volts. 

This estimate, however, neglects the possibility of a 
transient breakdown of the RC circuit approximation, 
prior to diffusive relaxation across the entire cell. In the 
case of response to a suddenly applied DC voltage, the 
salt adsorption by each the diffuse layer occurs over a 
time scale, Tc = AdL/Z?, during which diffusion spreads 
the nearby region of depleted neutral bulk solution over 
a distance 



'Dtc — y/XoL. Therefore, the requirement 
that the local bulk conductivity does not change signifi- 
cantly during charging dynamics yields the refined esti- 
mate m, 



ad 



cqvXdL 



< 1 



which replaces L/Xjj by y^L/X in the estimates above. 
This sets a lower bound for the limiting voltage for the 
validity of circuit models. 

In the case of AC forcing at frequency uj, the bound 
Q!c; ^ 1 is appropriate for low frequencies, 2itu!Tc <^ 1, 
but circuit models remain valid up to higher voltages at 
higher frequencies. At moderate frequencies, 27rwrc > 1, 
the double layer does not have enough time for complete 
charging, and the Warburg-like diffusion layer due to salt 
adsorption (which oscillates at twice the frequency) only 
propagates to a distance, y^D/ATruj. Therefore, we may 
crudely estimate {w / cq) Attuj / D ^ 1 to avoid signifi- 
cant changes in bulk conductivity in the diffusion layer. 



(A more careful estimate would take into account that 
only partial charging and salt adsorption occur with in- 
creasing frequency.) At higher frequencies, tuXjj/D « 1, 
the diffuse layer does not have enough time to equilibrate, 
and little charging occurs. 

These arguments can be made more precise using 
matched asymptotic expansions to describe the thin dou- 
ble layer limit, starting from an explicitly time-dependent 
model. For the case of one-dimensional response to a 
suddenly applied voltage, this was done in Ref. [6] start- 
ing from the Poisson-Nernst-Planck equations of (time- 
dependent) dilute solution theory. In part II, we will 
derive modified PNP equations for dynamics with steric 
effects and repeat the same kind of asymptotic analysis 
to reach a similar conclusion our simple arguments here: 
Steric effects greatly extend the range of applicability of 
RC circuit models, compared to what would be expected 
on the basis of dilute solution theory. 



C. Diffuse-layer surface conduction 



Another feature not present in circuit models is the 
possibility of current passed along the surface through the 
diffuse layer, demonstrated by Bikerman [T3| and consid- 
ered extensively in theories of electrokinetics by Dukhin 
and collaborato rs l3Cl| . As first shown in Ref. [61 and elab- 
orated in Ref. [2J] in the setting of dilute solution the- 
ory, the relative strengths of tangential "surface fluxes" 
through the diffuse layer, compared to bulk fluxes, are 
controlled by the same dimensionless groups that govern 
ion adsorption (discussed above). This is actually quite 
a general result, as we now briefly explain. We will thus 
conclude that steric effects also greatly reduce the impor- 
tance of surface conduction in the diffuse layer compared 
to the classical predictions of dilute solution theory. 

Assuming small local perturbations from thermal equi- 
librium, the flux density (number/area-time) of ionic 
species i is given by 



Fi = -biCi\7fj.i 



(37) 



where the chemical potential fj,i generally has a differ- 
ent form than ([1]) in a concentrated solution (e.g. see 
Part II). Consider a thin diffuse layer near a charged 
surface, where the ion concentration Ci departs from its 
nearby neutral bulk value c^. Due to fast relaxation at 
the small scale of the screening length, the diffuse-layer 
concentration remains in quasi-equilibrium at nearly con- 
stant chemical potential in the normal direction, /i^ ~ /i^, 
in the thin double layer limit. There can, however, be 
small tangential gradients, V||/Lti 7^ at the macroscropic 
length scale leading to an excess diffuse-layer "surface" 
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flux density (number/length-time): 

{-b^C^Wl\^I^+hc'',W\l^i^,)dx (38) 

-/ surface 

(39) 



(^bulk 
rfa 

pbulk 

-V||A*- / h{c,-c\)dx 

J surf ace 



For a constant mobility bi, this takes the same form as 
the bulk flux density ((37)) . 



(40) 



where the bulk concentration (number / volume) has been 
replaced by the diffuse-layer "surface concentration" 
(number /area) , 



bulk 



[a - c- ) dx 



(41) 



surface 



In a concentrated solution, we generally expect that the 
mobility bi might decrease in the difuse layer, due to 
steric effects and large normal electric fields, so this for- 
mula may overestimate the surface flux density. 

Following Bikerman [l^, [HI, [HI , we may estimate the 
relative importance of surface to bulk flux densities at a 
length scale L by the dimensionless group 



IL 

F^L 



(42) 



which we see also measures the relative importance of 
"surface adsorption" of ion i in the diffuse layer relative to 
the bulk concentration. For a highly charged diffuse layer 
I^dI 3> the ions are mostly of one type (counter- 
ions), so Fi ~ 5 ~ w. The surface current is also carried 
mostly by those ions, J* ~ zeFf, while the bulk current is 
carried by both ions, « 2{ze)Fj' (neglecting diffusion 
compared to electromigration). Therefore, we see that 
the "Dukhin number" comparing surface conduction at a 
highly charged diffuse layer to bulk conduction at nearly 
uniform concentration, 



Du 



J' 



JbL 2coL 



2 



(43) 



is roughly half of the dimensionless group as governing 
salt adsorption by the diffuse layer. 

We have seen that steric effects greatly reduce w com- 
pared to the predictions of dilute solution theory, so that 
as — 0{\d/L) remains small up to rather high volt- 
ages. Since the calculation above over-estimates the im- 
portance of surface conduction, it is clear that steric ef- 
fects also greatly reduce the Dukhin number compared 
to the predictions of dilute solution theory. We conclude 
that surface conduction in a thin diffuse layer does not 
become important until voltages large enough to violate 
the equivalent circuit approximation are applied across 
it. Compact layer surface conduction may still be im- 
portant in some cases, but it too is limited by the same 
steric effects. 



V. COMPACT LAYER EFFECTS 

We now check the robustness of these conclusions to 
the additional presence of an insulating surface layer be- 
tween the metal electrode and the electrolyte where the 
EDL models are applied. We suppose that this layer is 
not involved in the dynamics so that its properties do not 
change, and assume for the sake of simplicity that these 
properties are linear so that it can be described as a fixed 
surface capacitance Cg- 

Such models have been used in many circumstances, 
sometimes to describe the Stern layer corresponding to 
condensed ions. Our approach here is different in the 
sense that a layer of condensed ions would be involved in 
the charging/discharging process. What we have in mind 
is closer to a thin film, e.g. of oxide on the metal, of thick- 
ness hs and dielectric constant ts so that Cs = Cs/hs- 
This form has been proposed to model coating layers on 
electrodes in the context of AC electrokinetics [ij . 

The overall differential capacitance of the interface is 



1 

C 



1 



1 
'C~s 



with Cs — dq/d'^Sj where q is the total charge per area in 
the double layer, and ^'s is the voltage drop across the 
aforementioned compact layer. The total voltage drop 
across the interface is ^ = 'ifs + "^d, and is the 
total capacitance of the interface. 

A useful dimensionless parameter to quantify the ef- 
fect of this surface layer is ^ = e/{XDCs), which has 
been employed recently in general studies of diffuse- 
charge dynamics [a. |l 6| . as well as in theoretical mod- 
[3, [lo, 58^ 7a [94I and in the fitting of exper- 
l3, [731 for induced-charge (AC or DC) 
No surface layer corresponds to S = 0. 



els m, 

imental data 
electrokinetics. 
In the above mentioned picture of the layer of oxide, 
d — {e/ts){hs/XD) so that even very thin oxide layers 
can yield not too small values of ^ if e/cs is large. 

Let us now compare the differential capacitance C 
(with surface layer) to the bare double layer differential 
capacitance previously plotted in FiglTl for the PB equa- 
tion and the two models with steric effects. The corre- 
sponding plots are provided on FigdOlfor 5 = 0.25. 

Many qualitative points are obvious. First, the PB 
differential capacitance does not blow up exponentially 
anymore, as the surface layer takes over when the double 
layer voltage gets large, so that C converges to the finite 
value Cs within the PB model. Most of the potential 
drop is then across the surface layer. Second, the two 
main consequences of steric effects pointed above remain 
valid: C is weaker when steric effects are taken into ac- 
count, and C(4') is nonmonotonous with a maximum at 
intermediate values and a further decrease to zero. 

When the parameter 6 is zero or small(i.e. the Stern 
layer capacitance is large), a closely packed "diffuse" 
layer has to form in response to a high potential drop. 
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FIG. 10: The total double layer differential capacitance in 
presence of a Stern layer with S — 0.25 versus the total po- 
tential drop across the double layer. 
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FIG. 11: The differential capacitance C of the double layer as 
a function of potential difference 'i' for various values of Stern 
layer capacitance values. 



This is the case when the steric effects become impor- 
tant and the PB theory fails. However, for the higher 
values of (5, most of the potential drop is reahzed over 
the Stern layer, and the ions in the diffuse layer are not 
as densely packed. Thus the system stays below the steric 
limit, and the PB theory agrees with the MPB theory, as 
demonstrated by the total capacitance plots of Fig. [TT] 
Similar comments can be made on Fig llOi which shows 
that for 6 = 0.25, steric effects may or may not be impor- 
tant depending on value of the dimensionless parameter 
v = 2a^co. For the distribution of the voltage drop over 
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FIG. 12: The Stern layer voltage drop as a function of the 
total double layer voltage difference. For reference, the line 
^5 = is also drawn. In the PB theory, almost all of the 
voltage drop is realized over the Stern layer when the voltage 
drop is on the order of a few voltages. 



the Stern layer and the diffuse layers, see FigfT^ 

This analysis, although correct, is somewhat mis- 
leading, however, since the assumption of a constant 
compact-layer capacitance may not always be reason- 
able. In many situations of interest, where several volts 
(« lOOfcT/e) are applied across the double layer, it is 
unlikely that the compact layer could withstand a signif- 
icant fraction of the total voltage. Note that dielectric 
breakdown in water can occur in average fields as low as 
20 MV/m = 0.02 V/nm in experiments applying submi- 
crosceond voltage pulses (well below the charging time of 
the double layers) [5l|. The critical field may be higher 
in the Stern layer, where water is confined by ions at the 
outer Helmholtz plane fisll, but it seems implausible for 
an atomic layer to sustain several volts without some- 
how "short circuiting" via electron cascades, electrolysis, 
Faradaic reactions, etc. In some cases, the compact layer 
models a dielectric thin film, which may be considerably 
thicker. Again, however, most coating materials, such as 
teflon or various metal oxides, undergo dielectric break- 
down in fields of order 10 MV/m, so a dielectric coating 
cannot easily withstand several volts unless it is at least 
10 nm wide. 
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In general, we see that the capacitance of the compact 
layer must effectively decrease at large voltages, which 
corresponds to the limit 5^0. As a result, a signficant 
fraction of a large voltage must be sustained by the dif- 
fuse layer, making steric effects important in many situa- 
tions of interest. Regardless of the accuracy of our simple 
models, therefore, we believe that the predicted qualita- 
tive effects of ion crowding are likely to have broad rele- 
vance for experimental systems applying large voltages. 

VI. CONCLUSION 

We have used two simple models for the double layer to 
account for crowding effects which necessarily take place 
at intermediate and large applied voltages. These mod- 
els are both based on modifications of the Poisson Boltz- 
mann description of dilute solutions. This strategy has 
lead us to identify important operational consequences of 
these crowding effects the thin double layers, namely a 
largely reduced double layer capacitance and a decreased 
ion uptake from the bulk. 

We have provided in these sections explicit formulas 
for the total and differential capacitance and for the salt 
uptake of interfaces at equilibrium, as a function of the 
potential drop across the interface. More precisely we 
have recalled the PB results (no steric effects v = 0), and 
given results for the two simple models with steric effects 
{ly ^ 0), considering for all cases the possibility of a finite 
insulating layer on the electrode {S ^ 0). 

The two models lead to remarkably similar results sug- 
gesting that these semi-quantitatively hold beyond the 
specifics of these models. Both showed marked differ- 



ences with the PB approach: the differential capacitance 
and salt uptake are much weaker, and the former varies 
non-monotonously with the applied potential. 

These observations for the equilibrium properties have 
led us to make predictions for the dynamics of electrolyte 
cells of size quite larger than the Debye length: (i) an ef- 
fective equivalent RC circuit description holds for a wider 
range of potentials than expected on the simple basis of 
the PB equation, (ii) the response time is much smaller 
than expected from PB at large voltages, (iii) this time 
decreases at large voltages after an initial increase for 
lower values. 

The dramatic effect of steric constraints in this prob- 
lem also shows that other predictions of nonlinear PB 
theory, such as the change of scaling from to 
\V\ log \ V\ for AC electro-osmosis [111, are limited in ap- 
plicability and should be revisited with models taking 
crowding effects into account. More generally, this work 
suggests that, beyond the present problem of ionic trans- 
port in electrolytic systems, the description of electroki- 
netic effects at large applied voltages should be revisited 
to correct shortcomings of dilute-solution theory. 

Acknowledgments 

This work was supported in part by the MRSEC pro- 
gram of the National Science Foundation under award 
number DMR 02-12383 (MZB). A A acknowledges very 
gratefully the hospitality of the Department of Mathe- 
matics at MIT where this study was realized. 



[1] A. Ajdari. AC pumping of liquids. Phys. Rev. E, 61:R45- 
R48, 2000. 

[2] S. Alexander, P. M. Chaikin, P. Grant, G. J. Morales, 
and P. Pincus. Charge renormalization osmotic pressure 
and bulk modulus of colloidal crystals: theory. J. Chem. 
Phys., 80:5776-5781, 1984. 

[3] M. Z. Bazant and Y. Ben. Theoretical prediction of fast 
3d ac electro-osmotic pumps. Lab on a Chip, 6:1455- 
1461, 2006. 

[4] M. Z. Bazant, K. T. Chu, and B. J. Bayly. Current- 
voltage relations for electrochemical thin films. SIAM J. 
Appl. Math., 65:1463-1484, 2005. 

[5] M. Z. Bazant and T. M. Squires. Induced-charge electro- 
kinetic phenomena: Theory and microfluidic applica- 
tions. Phys. Rev. Lett, 92:066101, 2004. 

[6] M. Z. Bazant, K. Thornton, and A. Ajdari. Difi'use 
charge dynamics in electrochemical systems. Phys. Rev. 
E, 70:021506, 2004. 

[7] Martin Z. Bazant and Yuxing Ben. Theoretical prediction 
of fast 3d ac electro-osmotic pumps. Lab on a Chip, 2006. 

[8] L. Belloni. A hypernetted chain study of highly asym- 
metrical poly-electrolytes. Chem. Phys., 99:43, 1985. 

[9] L. Belloni. Ionic condensation and charge renormaliza- 



tion in colloidal suspensions. Colloids and Surfaces A: 
Physicochemical and Engineering Aspects, 140:227-243, 
1998. 

[10] J. J. Bikerman. Z. Phys. Chem. Abt. A, 163:378, 1933. 

[11] J. J. Bikerman. Electrokinetic equations and surface 
conductance, a survey of the difltuse double layer theory 
of colloidal solutions. Trans. Faraday Soc, 36:154-160, 
1940. 

[12] Victor A. Bloomfield. Dna condensation by multivalent 
cations. Biopolymers, 44:269-282, 1998. 

[13] J. O.'M Bockris and A. K. N. Reddy. Modern Electro- 
chemistry. Plenum, New York, 1970. 

[14] K. Bohinc, A. Iglic, T. Slivnik, and V. Kralj-Iglic. 
Charged cylindrical surfaces: EflFect of finite ion size. Bio- 
electrochemistry, 57:73-81, 2002. 

[15] K. Bohinc, V. Kralj-Iglic, and A. Iglic. Thickness of elec- 
trical double layer, effect of ion size. Electrochim. Acta, 
46:3033-3040, 2001. 

[16] A. Bonnefont, F. Argoul, and M.Z. Bazant. Analysis of 
difi'use-layer efi^ects on time-dependent interfacial kinet- 
ics. J. Electroanal. Chem., 500:52-61, 2001. 

[17] I. Borukhov. Charge renormalization of cylinders and 
spheres: Ion size effects. J. Pol. Set. B: Pol. Phys., 



16 



42:3598-3620, 2004. 
[18] 1. Borukhov, D. Andclman, and H. Orland. Steric effects 

in electrolytes: A modified poisson-boltzmann approach. 

Phys. Rev. Lett, 79:435-438, 1997. 
[19] I. Borukhov, D. Andelman, and H. Orland. Adsorption 

of large ions from am electrolyte solution: A modified 

poisson boltzmann equation. Electrochim. Acta, 46:221- 

229, 2000. 

[20] A. B. D. Brown, C. G. Smith, and A. R. Rormic. Pump- 
ing of water with ac electric fields applied to asymmetric 
pairs of microelectrodes. Phys. Rev. E, 63:016305, 2000. 

[21] S. L. Carnie and G. M. Torrie. The statistical mechanics 
of the electrical double- layer. Adv. Ghem. Phys., 46:141- 
253, 1984. 

[22] D. L. Chapman. A contribution to the theory of electro- 
capillarity. Philos. Mag., 25:475-481, 1913. 

[23] K. T. Chu and M. Z. Bazant. Electrochemical thin films 
at and above the classical limiting current. SIAM J. Appl. 
Math., 65:1485-1505, 2005. 

[24] K. T. Chu and M. Z. Bazant. Nonlinear electrochem- 
ical relaxation around conductors. Physical Review E, 
74:011501, 2006. 

[25] E. Darve D. Saintillan and E. Shaqfeh. Journal of Fluid 
Mechanics, 2006. 

[26] P. Debye and H. Falkenhagen. Physik. Z., 29:121, 1928. 

[27] B. V. Deryagin and S. S. Dukhin. The theory of surface 
conductance. Colloid. J. USSR., 31:277-283, 1969. 

[28] N. J. Dudnoy, J. B. Bates, D. Lubben, and F. X. Hart. 
Thin-film rechargeable lithium batteries with amorphous 
Lia;Mn204 cathodes. In J. Bates, editor. Thin Film Solid 
Ionic Devices and Materials, pages 201-214. The Elec- 
trochemical Society, Pennington, NJ, 1995. 

[29] S. S. Dukhin. Electrokinetic phenomena of the second 
kind and their applications. Adv. Colloid Interface Sci., 
35:173-196, 1991. 

[30] S. S. Dukhin. Non-equilibrium electric surface phenom- 
ena. Adv. Colloid Interface Sci., 44:1-134, 1993. 

[31] M. Eigen and E. Wicke. The thermodynamics of elec- 
trolytes at higher concentration. J. Phys. Chem., 58:702- 
714, 1954. 

[32] V. Freise. Zur theorie der diffusen doppelschicht. Z. 
Elektrochem., 56:822-827, 1952. 

[33] J.B. Freund. Electroosmosis in a nano-meter scale chan- 
nel studied by atomistic simulation. J. Chem. Phys., 
116:2194 2200, 2002. 

[34] A. Frunikin. Wasserstoffiiberspannung und struktur der 
doppelschict. Z. Phys. Chem., 164A:121-133, 1933. 

[35] N. I. Gamayunov, V. A. Murtsovkin, and A. S. Dukhin. 
Pair interaction of particles in electric field. 1. features of 
hydrodynamic interaction of polarized particles. Colloid 
J. USSR, 48:197-203, 1986. 

[36] L. A. Geddes. Historical evolution of circuit models for 
the electrode-electrolyte interface. Ann. Biomedical Eng., 
25:1-14, 1997. 

[37] D. Gillespie, W. Nonner, and R. S. Eisenberg. Coupling 
poisson-nernst-planck and density functional theory to 
calculate ion fiux. J. Phys.: Condens. Matter, 14:12129- 
12145, 2002. 

[38] D. Gillespie, W. Nonner, and R. S. Eisenberg. Density 

functional theory of charged, hard-sphere fluids. Phys. 

Rev. E, 68:031503, 2003. 
[39] D. Gillespie, M. Valisko, and D. Boda. Density functional 

theory of the electrical double layer: the rfd functional. 

J. Phys.: Condens. Matter, 17:6609-6626, 2005. 



[40] A. Gonzalez, A. Ramos, N. G. Green, A. Castellanos, and 
H. Morgan. Fluid flow induced by non-uniform ac elec- 
tric flelds in electrolytes on microelectrodes. ii. a linear 
double-layer analysis. Phys. Rev. E, 61:4019, 2000. 

[41] M. Gouy. Sur la constitution de la charge electrique a la 
surface d'un electrolyte. J. de Phys., 9:457-468, 1910. 

[42] L. Guldbrand, B. Jonsson, H. Wennerstrom, and P. J. 
Linse. Electrical double-layer forces: A monte carlo 
study. J. Chem. Phys., 80:2221, 1984. 

[43] B.-Y. Ha and Andrea J. Liu. Counteriou-niodiatcd at- 
traction between two like-charged rods. Phys. Rev. Lett., 
79:1289-1292, 1997. 

[44] H. Helmholtz. Studien fiber electrische grenzschichten. 
Ann. Phys. Chem., 7 (ser. 3):337-382, 1879. 

[45] J. R. Henderson, editor. Fundamentals of Inhomogeneous 
Fluids. Marcel Dekker, New York, 1992. 

[46] M. L. Henle, C. D. Santangelo, D. M. Patel, and 
P. A. Pincus. Distribution of counterions near discretely 
charged planes and rods. Europhys. Lett, 66:284-290, 
2004. 

[47] R. J. Hunter. Foundations of Colloid Science. Oxford 
University Press, Oxford, 2001. 

[48] A. Iglic and V. Kralj-Iglic. Influence of finite size of ions 
on electrostatic properties of electric double layer. Elec- 
trotecnical Rev. (Slovenia), 61:127-133, 1994. 

[49] D. Henderson J. Reszko-Zygmunt, S. Sokolowski and 
D. Boda. Temperature dependence of the double layer 
capacitance for the restricted primitive model of an elec- 
trolyte solution from a density functional approach. J. 
Chem. Phys., 122:084504, 2005. 

[50] L. Joly, C. Ybert, E. Trizac, and L. Bocquet. Hydrody- 
namics within the electric double layer on slipping sur- 
faces. Phys. Rev. Lett, 93:257805, 2004. 

[51] H. M. .Jones and E. E. Kunliardt. Pulsed dielectric break- 
down of pressurized water and salt solutions. Journal of 
Applied Physics, 77:795-805, 1995. 

[52] Bo Jonsson, A. Nonat, C. Labbez, B. Cabane, and 
H. Wennerstrom. Controlling the cohesion of cement 
paste. Langmuir, 21:9211-9221, 2005. 

[53] D. Kim and E. Darve. Molecular dynamics simulations 
of rough wall microchannels. Phys. Rev. E, 73:051203, 
2006. 

[54] V. Kralj-Iglic and A. Iglic. A simple statistical mechan- 
ical approach to the free energy of the electric double 
layer including the excluded volume effect. J. Phys. 11 
France, 6:477-491, 1996. 

[55] C.W. Outhwaite L.B. Bhuiyan and D. Henderson. A 
modified poisson-boltzmann analysis of the capacitance 
behaviour of the electric double layer at low tempera- 
tures. J. Chem. Phys., 123:034704, 2005. 

[56] S. Levine and C.W. Outhwaite. Comparison of theories 
of the aqueous electric double layer at a charged plane in- 
terface. J. Chem. Sac, Faraday Trans. 2, 74:1670, 1977. 

[57] J. A. Levitan, S. Devasenathipathy, V. Studer, Y. Ben, 
T. Thorsen, T. M. Squires, and M. Z. Bazant. Ex- 
perimental observation of induced-charge electro-osmosis 
around a metal wire in a microchannel. Colloids and Sur- 
faces A, 267:122-132, 2005. 

[58] H. Bruus L.H. Olesen and A. Ajdari. Ac electroki- 
netic pumps: The effect of geometrical confinement, 
faradaic current injection and nonlinear surface capac- 
itance. Phys. Rev. E, 73:056313, 2006. 

[59] J. Lyklema. Fundamentals of Interface and Colloid Sci- 
ence. Volume I: Fundamentals. Academic Press Limited, 



17 



San Diego, CA, 1991. [79 
J. Lyklema. Fundamentals of Interface and Colloid Sci- 
ence. Volume II: Solid-Liquid Interfaces. Academic Press 
Limited, San Diego, CA, 1995. [80 
J. Lyklema. Does electrical double layer formation lead 
to salt exclusion or to uptake? Phys. Rev. E, 71:032501, [81 
2005. 

J. R. Macdonald. Theory of the differential capacitance [82 
of the double layer in unadsorbed electrolytes. J. Chem. 
Phys., 22:1857-1866, 1954. 

J. R. Macdonald. Impedance spectroscopy: Old problems [83 
and new developments. Electrochim. Acta, 35:1483-1492, 
1990. 

G. S. Manning. Limiting laws and counterion condensa- [84 
tion in polyelectrolyte solutions i. colligative properties. 

J. Chem. Phys., 51:924-933, 1969. 

S. Marcelja. Electrostatics of membrane adhesion. Bio- 
phys. J., 61:1117, 1992. [85 
V. A. Murtsovkin. Nonlinear flows near polarized dis- 
perse particles. Colloid Journal, 58:341-349, 1996. [ 
B. Nadler, Z. Schuss, A. Singer, and R.S. Eisenberg. 
Diffusion through protein channels: from molecular de- 
scription to continuum equations. Nanotechnology, 3:439, [87 
2003. 

B. Nadler, Z. Schuss, A. Singer, and R.S. Eisenberg. Ionic 
diffusion through confined geometries: from langevin [. 
equations to partial differential equations. J. Phys.: Con- 
dens. Matter, 16:S2153-S2165, 2004. 

B. J. Neudecker, N. J. Dudney, and J. B. Bates, "lithium- [89 
free" thin-film battery with in situ plated Li anode. J. 
Electrochem. Soc, 147:517-523, 2000. 

J. Newman. Electrochemical Systems. Prentice-Hall, Inc., 
Englewood Cliffs, NJ, second edition, 1991. 

F. Oosawa. Interaction between parallel rodlike [90 

macroions. Biopolymers, 6:134, 1968. 

F. Oosawa. Polyelectrolytes. Marcel Dekker, New York, 

1971. [91 

C. W. Outhwaite and L.B. Bhuiyan. Theory of the electric 
double layer using a modified poisson-boltzmann equa- 
tion. J. Chem. Soc, Faraday Trans. 2, 76:1388-1408, [92 
1980. 

C.W. Outhwaite and L.B. Bhuiyan. A further treatment [93 
of the exclusion-volume term in the modified poisson- 
boltzmann theory of the electric double layer. J. Chem. 
Soc, Faraday Trans. 2, 78:775-785, 1982. [94 
C.W. Outhwaite and L.B. Bhuiyan. An improved mod- 
ified poisson-boltzmann equation in electric-double-layer [95 
theory. J. Chem. Soc, Faraday Trans. 2, 78:707-718, 
1983. 

R. Qiao and N. R. Aluru. Ion concentrations and velocity [96 
profiles in nanochannel electroosmotic flows. J. Chem. 
Phys., 118:4692-4701, 2003. 

A. Ramos, , H. Morgan, N. G. Green, A. Gonzalez, and 

A. Castellanos. Pumping of liquids with traveling-wave [97 

electro-osmosis. Journal of Applied Physics, 97:084906, 

2005. 

A. Ramos, A. Gonzalez, A. Castellanos, N. G. Green, and 

H. Morgan. Pumping of liquids with ac voltages applied [' 
to asymmetric pairs of microelectrodes. Phys. Rev. E, 
67:056302, 2003. 



A. Ramos, H. Morgan, N. G. Green, and A. Castellanos. 
AC electric-field-induced fluid flow in microelectrodes. J. 
Colloid Interface Set., 217:420-422, 1999. 

R. Roth and D. Gillespie. Physics of size selectivity. Phys. 
Rev. Lett., 95:247801, 2005. 

I. Rubinstein. Electro- Diffusion of Ions. SIAM Studies 
in Applied Mathematics, SIAM, Philadelphia, PA, 1990. 
I. Rubinstein and L. Shtilman. Voltage against current 
curves of cation exchange membranes. J. Chem. Soc 
Faraday. Trans. II, 75:231-246, 1979. 
I. Rubinstein and B. Zaltzman. Electro-osmotically in- 
duced convection at a permselective membrane. Phys. 
Rev. E, 62:2238-2251, 2000. 

D. J. MitcheU S. L. Carnie, D.Y.C. Chan and B.W. Nin- 
ham. The structure of electrolytes at charged surfaces: 
The primitive model. J. Chem. Phys., 74:1472-1478, 
1981. 

Christian D. Santangelo. Computing counterion densities 
at intermediate coupling, cond-mat/0501498 
Z. Schuss, B. Nadler, and R.S. Eisenberg. Derivation 
of pnp equations in bath and channel from a molecular 
model. Physical Review E, 64:036116, 2001. 
V. N. Shilov and T. S. Simonova. Polarization of electric 
double layer of disperse particles and dipolophoresis in a 
steady (dc) field. Colloid J. USSR, 43:90-96, 1981. 
I. N. Simonov and V. N. Shilov. Theory of low-frequency 
dielectric-dispersion of a suspension of ideally polarizable 
spherical particles. Colloid J. USSR, 39:775-780, 1977. 
T. A. Simonova, V. N. Shilov, and O. A. Shramko. 
Low-frequency dielectrophoresis and teh polarization in- 
teraction of uncharged spherical particles with an in- 
duced deby atmosphere of arbitrary thickness. Colloid 
J., 63:108-115, 2001. 

M. Sluyters-Rehbach and J. H. Sluyters. Electroanalytical 
Chemistry, volume 4, pages 1-128. Marcel Dekker, New 
York, 1970. 

W. H. Smyrl and J. Newman. Double layer structure at 
the limiting current. Trans. Faraday Soc, pages 207-216, 
1967. 

T. M. Squires and M. Z. Bazant. Induced-charge electro- 
osmosis. J. Fluid Mech., 509:217-252, 2004. 
T. M. Squires and M. Z. Bazant. Breaking symmetries 
in induced-charge electro-osmosis and electrophoresis. J. 
Fluid Mech., 560:65-101, 2006. 

O. Stern. Zur theorie der electrolytischen doppelschicht. 
Z. Elektrochem., 30:508-516, 1924. 

N. Takami, T. Ohsaki, H. Hasabe, and M. Yamamoto. 
Laminated thin Li-ion batteries using a liquid electrolyte. 
J. Electrochem. Soc, 148:A9-A12, 2002. 

B. Wang, J. B. Bates, F. X. Hart, B. C. Sales, R. A. 
Zuhr, and J. D. Robertson. Characterization of thin-film 
rechargeable lithium batteries with lithium cobalt oxide 
cathodes. J. Electrochem. Soc, 143:3204-3213, 1996. 

E. Wicke and M. Eigen. ber den einflu des raumbe- 
darfs von ionen in wssriger Isung auf ihre verteilung im 
elektrischen feld und ihre aktivittskoeffizienten. Z. Elek- 
trochem., 56:551, 1952. 

E. Yariv. Induced-charge electrophoresis of nonspherical 
particles. Phys. Fluids, 17:051702, 2005. 



